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We studied complex spectra of spin-two boson systems represented by E®e and E®(&i + 62) Jahn- Teller models. 
For E®e, at particular rotation quantum numbers we found a coexistence of up to three regions of the spectra, 
(i) the dimerized region of long-range ordered (extended) pairs of oscillating levels, (ii) the short-range ordered 
(localized) "kink lattice" of avoiding levels, and (iii) the intermediate region of kink nucleation with variable 
range of ordering. This structure appears above certain critical line as a function of interaction strength. The 
level clustering and level avoiding generic patterns reflect themselves in several intermittent regions between up-to 
three branches of spectral entropies. Linear scaling behavior of the widths of curvature probability distributions 
provides the conventionally adopted indication for the presence of quantum chaos. The mapping onto classical 
integrable Calogero-Moser gas provided useful insight into the complex level dynamics, including the soliton 
collisions representing the level avoidings and, in a range of model parameters, a novel view on the notion of 
quantum chaos formulated in terms of quantum numbers via the logistic equation. We found that apart from two 
limiting cases of J3<8>(6i + 62) model {E®e and Holstein model) the distribution of nearest neighbor spacings of this 
model is rather stable as to the change of parameters and different from Wigner one. This limiting distribution 
assumably shows scaling ~ vS at small 5 and resembles the semi-Poisson law P(S) = 45exp( — 25) at 5 > 1. 
The latter is believed to be universal and characteristic, e.g. at the transition between metal and insulator phases. 



1. Introduction 

The results of intensive research in quantum 
mechanics show (^-[3]) that the smooth wave- 
like quantum world contains elements of quan- 
tum complexity and that symptoms of chaos 
enter even into the wave patterns associated 
with atomic energy levels. Starting with nu- 
clear physics 0, one observes complexity in 
quantum phenomena in spectral features, eigen- 
functions, dynamics of molecular processes, solid 
state physics and information science (Berry 0] 
suggested the "quantum chaology" being the 
adequate term for a quantum system whose 
classical counterpart exhibits chaos instead the 
mostly used and widely accepted term " quantum 
chaos"). 

Besides the quantum systems that exhibit chaos 



in the classical or a semiclassical limit there ex- 
ists a class of quantum systems that do not have 
a naive classical or semiclassical counterpart. In 
addition, the way of passing to semiclassical ap- 
proximation, e.g., in spin-boson systems is not 
unique and presents essential ambiguity due to 
different possible ways of decoupling quantum 
variables [2]. Different ways of performing semi- 
classical approximation are known to lead to dif- 
ferent answers concerning the chaotic behaviour 
of the system. This ambiguity means that a clas- 
sical analogue of such a model is not well-defined 
and cannot be a reliable object in exploring the 
quantum chaos issues. 

The class of "spin-boson" models, e.g., two and 
more levels coupled to bosons (vibrons, phonons, 
photons) were reported to exhibit chaotic phe- 
nomena first by Lewenkopf et al [Hj and Cibils et 
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al and many others. Fujisaki et al (| 9I1U| ) in- 
vestigated the systems with the absence of a rea- 
sonable quasiclassical limit. There such phenom- 
ena as quantum tunneling and multi-mode nona- 
diabatic fluctuations lead to complicated wave 
packet dynamics which breaks down the Born- 
Oppenheimer adiabatic approximation. 
Typical representative of such a non-adiabatic 
system is the two-level molecular system cou- 
pled with two vibron modes of different symme- 
try against the transformation of reflection - the 
Jahn- Teller (JT) class of models. The chaos in JT 
molecules was investigated within a generalized 
model by Cederbaum et al in a series of papers 
Recently, Yamasaki et al f° r the nrs t 
time investigated a possibility of chaos in spectra 
of the E(g>e JT model. Their analysis was based 
on the approximation of the Hamiltonian by sep- 
arating it into the adiabatic "Mexican hat" po- 
tential and an additional part supplemented via 
a parameter for a nonlinear mode coupling of a 
trigonal symmetry to include the effect of fluctua- 
tions and non-integrability. Thus, the nonlinear- 
ity was included via mode coupling in addition 
to the mean field bare part of the Hamiltonian. 
The authors concluded that the quantum chaos 
reflected itself in the Wigner level spacing distri- 
bution as a consequence of the said nonlinearity 
of the Hamiltonian, meanwhile for the linear part 
the absence of these patterns was stated. 
In the light of said problems due to the absence 
of a reasonable quasiclassical approximation we 
have used an exact numerical approach to the 
characteristics of excited (quasi-continuum) spec- 
tra of the E®e JT model |13I14| . This approach 
differs from that by Yamasaki et al. in the fol- 
lowing: a) we do not introduce an explicit non- 
linearity into the initial Hamiltonian. However, 
its SU(2) symmetry involves an intrinsic nonlin- 
earity which is revealed by the exact elimination 
of the electronic degrees of freedom ^5]. b) We 
started from the E®(£>i + 62) model with differ- 
ent coupling strengths for both modes. The rota- 
tion symmetric E(8>e model represents its partic- 
ular case with equal coupling constants and thus 
with the symmetry of higher (rotational) sym- 
metry group [the difference of the coupling con- 
stants in realistic systems is likely to be caused, 



for example, by the spatial anisotropy of crys- 
tals] . The importance of starting from the more 
general symmetry is evident from the nature of 
quantum fluctuations: namely, the variational ap- 
proach to E(g>(6i + 62) model used for the cal- 
culation of the ground state \H>\ was shown to 
yield the largest deviations just for the rotation 
symmetric case. In the ground state an abrupt 
change of the energy at equal couplings is an 
artefact of the adiabatic approximation: the en- 
ergy region of the quantum E(g>e model is situ- 
ated within a smeared border around /3/a = 1 
(a and (3 being the coupling strengths of the an- 
tisymmetric and the symmetric modes, respec- 
tively) between the " selftrapping" (a > (i) and 
" tunneling" (a < (3) part of the ground state of 
the model with broken rotational symmetry. In 
the smeared transition region both phases coex- 
ist and are quantum correlated (entangled) via 
nonlinear correlations between the modes. There 
the phonon-assisted (by the symmetric mode) 
tunneling contribution to the energy [from the 
admixture of two reflection symmetric levels in 
the excitation reflection Ansatz for the varia- 
tional wave function |15| ] results in the essential 
decrease of the ground state energy because of 
changing parity of the wave function by the par- 
ity reversal of the phonon 1-coordinates by the 
operator R, Rb x = -b x R, (3{^\Q 2 R\^\~) = 
(3^+\Q 2 \^+) » 0{*+\Q 2 \*i~). Therefore, it 
has its consequence in a dramatic increase of the 
Ham factor and in the corresponding decrease of 
the energy. 

Excited states follow the structure of dominat- 
ing selftrapping (a) or tunneling (/3) interactions. 
The mixing of the phases due to the tunneling is 
amplified by the analogous to the above matrix 
element via the reflection change of the parity in 
matrix elements with excited wave functions. Re- 
spective decrease of the energy manifests itself in 
the avoiding of the excited energy levels. 
The phase transition at a = f3 to the rotational 
symmetry phase provides a representation by the 
additional rotational quantum number j. At 
given j in E(g>e model the analogous phonon as- 
sisted tunneling between adjacent excited levels 
occurs and is responsible for the flip (kink) be- 
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tween the levels (this mechanism can be well rec- 
ognized, e.g., in Fig. lb). This nonlinear effect 
causes a mixture of the adjacent (even and odd) 
levels at stronger couplings and leads to the nucle- 
ation of the new kink lattice phase with increasing 
coupling constant. The tunneling mechanism due 
to short-range quantum fluctuations at large j is 
then the origin of peculiarly irregular behaviour 
of the related level spacing distributions and does 
not allow for full developing of the Wigner chaos. 
The phonon assisted tunneling and nucleation of 
the kink lattice phase can be also demonstrated 
within the formalism of a Calogero-Moser gas of 
pseudo-particles with repulsive interaction 



2. Hamiltonian approach: two-level-two- 
boson Hamiltonians 

The aim of this Section is to explore the nu- 
merical solution of the E<g>e Hamiltonian for the 
excited energy spectrum and wave functions in 
search for the statistical evidences of the symp- 
toms of chaotic behaviour at the quantum level. 
Namely, the participation ratios, spectral en- 
tropies and widths of level curvature distributions 
are anticipated to provide the relevant indications 
thereof. 

The E®e Hamiltonian with the parameter of non- 
integrability a is written as 

H = {b\b 1 +blb 2 +l)I+oi{b\+b 1 )a z ^a{bl+b 2 )a x , (1) 

where aj are Pauli matrices, / is the unit matrix. 
The pseudospin notation is used for the represen- 
tation of a two-level one-electron (spinless) sys- 
tem. 

For the secular equations of the vibronic parts 
of the gerade solution K = +1 (j — 1/2 = 
0, 2, 4, 6, ... ; —2, —4, . . .) the representation of ra- 
dial coordinates 



• ^ (.1 ■ id I)*) h + aV2rf k = Ef t 



will be used, and radial wave functions are as- 
sumed in the form f — c„$„, where 



2r 2 



r ~ 2rdrVdr + 2 r 



i^k, i,k= 1,2 (2) 
(3) 



$ + (n r , j, r), n = 2n r , 

<£ _ (n r ,j, r), n — 2n r + 1 . 



(4) 



The only non-zero off-diagonal elements are 
f nrnr = V2 J $+ $" r r ■ rdr = V2^n r + 1^ 

fn,, + ln r = V2 J ^tr+l^nJ ' rdr = -VV™r + 1 ■ 

The secular equations are written in the usual 
tridiagonal form (E^ is the "unperturbed" en- 
ergy for corresponding n,j), 

c n E^ + a(c„+i/„ 

n = 2n r , 

c n E^ +a(c n - 1 f n lfn r + ln r ) — Ec n , 

n = 2n r + 1 . (7) 

Numerical solutions of (7J| for the spectra and 
corresponding wave functions are exemplified in 
Figs. (1-4). 

3. Statistical methods: Level spacing and 
level curvature probability distribu- 
tions, long-range statistics 

The multitude of avoided crossings of energy 
levels is generally claimed to be a testimony of 
quantum chaos. The fundamental problem of 
the quantum-mechanical manifestation of classi- 
cal chaotic motion was investigated by statisti- 
cal methods based on the random matrix theory 
(RMT). It setforth reference patterns to which 
the presumably quantum chaotic systems should 
conform |2I16| . The conventionally adopted sta- 
tistical measures are either the short- or long- 
range statistics. The common representative of 
the former is the distribution of nearest neigh- 
bour level spacings (NNS), as well as the statis- 
tics involving the "curvatures" K = d 2 E n /da 2 
as functions of the nonintegrability parameter 
a. Among the latter a generally adopted tool 
is the A3-statistics [IS]. Both short- and long- 
range statistical measures have reference formu- 
las provided within the frames of the RMT. Thus, 
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Figure 1. Energy spectra for j = 19/2,255/2 as 
function of a. 

6x n a) a=l, j=255/2 S b ' a=1 - j=255/2 




Figure 2. Numerical wavefunctions \ n in the 
plane Qi x Q 2 , a = (3 = 2. States n=61, 63 
are localized (" exotic" ) . 



Figure 4. Reduced energies and spectral entropies 
of wavefunctions %n for different couplings a and 
j- 
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for the distribution of NNS the prediction for 
the fully developed chaos is the Wigner surmise 
P(S) ~ 5exp(— 5 2 ), for the opposite situation of 
the regular NNS distribution the Poisson distri- 
bution is representative P(S) = exp(— 5) (in the 
above formulas the spectra are supposed to be 
scaled so that (5) = 1 locally). The RMT leads to 
the indication of universality of level fluctuations 
conditioned solely by the time-reversal symmetry 
of the underlying Hamiltonians: the fluctuation 
patterns show a transition from the Poisson type 
to the Gaussian (orthogonal, unitary symplectic) 
type as the corresponding classical system shifts 
from the integrable to the chaotic regime. 

3.1. Spin one-boson model 

The class of non-integrable spin-boson models 
is known to exhibit quantum chaotic symptoms 
for one-boson many-level (2j + 1, j = 1/2,3/2,...) 
cases. Respective Hamiltonian 

H = ujb' t b + uj Q <j z + X[cr + b + cr^b t + 

e(a + tf + a_6)] , 
f± = (T x ±a y (8) 

is parametrized by the nonintegrability param- 
eter e. For e = 1 we obtain the exciton (Hol- 
stein) model (non-integrable, with evidence for 
chaos H]-|H]), while for e = one arrives at the 
Jaynes-Cummings model which is integrable. 

3.2. E(g)(6i+6 2 ) and E®e Jahn-Teller Model 

If a /3 but \a — f3\ <C a in (JIJ, there appear 
correlations between the excited levels with j =/= j 
of the Hamiltonian in radial representation due to 
the term ^2] 

<„-«,/*. • m 

Respective additional avoided level crossings 
change dramatically the probability distributions 
of level spacings when compared to those of the 
E ® e model (a = (3), Fig. 5. 

Fig. 5 presents the typical forms of the NNS 
distributions for different a, j3 of the JT model 
with broken radial symmetry. It is seen that for 
the special cases of higher symmetry the NNS 
distributions indeed are close to the Poisson dis- 
tribution (first row). Meanwhile far from the 



P(S) a=2 /3 = 2 P (S) a = 5 /3=1 




Figure 5. Sample nearest-neighbor distributions 
of levels (unfolded and scaled to (5) = 1) for 
different values of a, f3. The curves on the his- 
tograms represent Poisson (long dashed) , Wigner 
(full), semi-Poisson (stars) and y/S (crosses) dis- 
tributions. 



cases a — f3 and a <C (3, a ^> ft the NNS dis- 
tributions do not show a tendency to the RMT 
prediction (Wigner surmise). The well defined 
limit for these distributions seems not to have 
an analogue in the present versions of the RMT 
ensembles and presents some intermediate statis- 
tics about which we can make now only plausible 
assumptions. Thus, it is found that the semi- 
Poisson distribution P(S) = ASexp(-S) |17ll8] 
describes well this intermediate statistics in the 
domain of large S > 1 although for small S the 
limiting distribution seems to scale as P(S) ~ S u 
with v ranging between 0.3 -j- 0.5. The semi- 
Poisson statistics has recently found its applica- 
tion as a limiting case of level statistics at the M-I 
transition in Anderson localization models and, 
therefore, such an analogy seems to be a promis- 
ing item. As the second fitting formula (see 
Fig. 5, second row) we use the \/5-distribution 
P(S) ~ v^exp (— 3/25) which conforms to the 
interpolation of Brody type. From the numerical 
analysis of the nearest-neighbour level spacing 
probability distributions P(S) we conclude the 
diagnosis of the quantum chaotic patterns of JT 
models: 
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Figure 6. Standard deviations a 2 — (S* 2 ) — 1 of the 
NNS distributions. Grid lines correspond to the 
semi-Poisson (0.5) and \/~S (2/3) values. Far from 
the higher symmetry cases (Holstein and E®e JT 
model) the dispersion tends to a well defined limit 
a 2 ~ 0.5. The fully chaotic Wigner NNS distri- 
bution would have a 2 ~ 0.28. 



1. For small S, P{S) ~ S u , where v ac- 
quires a non-universal value from the inter- 
val (0.3,0.5). 

2. Except for small S, the limit cases of E ® 
(&i + 62) model of higher symmetry (E(g>e 
and Holstein) show up the close-to-Poisson 
distribution of NNS, especially in tails (Fig. 
5, first row). 

3. Far from the higher symmetry cases the 
NNS distribution tends to a universal law 
which is however markedly distinct from the 
Wigner distribution expected for the fully 
developed quantum chaos (Fig. 5, second 
row). 

4. The dispersion tends to a well defined and 
universal limit a 2 ~ 0.5 (Fig. 6) referring 
to the semi-Poisson distribution P(S) — 
4S'exp(— 2S). This distribution is known as 
a critical distribution at the M-I transition 
in Anderson localization models. 

For the E £g> e model a more pronounced in- 
dication for the presence of quantum chaos than 
the NNS level spacing distribution is the conven- 
tionally adopted |20j linear scaling behaviour of 
the widths of curvature probability distributions 




Figure 7. The widths of curvature distributions 
AK = (A 2 K) as function of j (a) and a (b). The 
scaling laws AK ~ j" and AK ~ a v are seen. 



AK = (A 2 K) as a function of nonintegrability 
parameter (a), Fig. 7. The quantum number j 
is a measure of a departure from the semiclassi- 
cal limit j — ► 00. From Fig. 8 it is seen that 
the widths of the curvature distribution show a 
scaling AK ~ j^a v with /1 ~ 1 and v ~ 0.5. 
The latter observation, as well as the form of the 
typical profiles of P(K) allows for the interpre- 
tation of the evolution of P(K) in the pseudo- 
time a as a solution of an ordinary diffusion-like 
equation with added ballistic ("telegraph") term 
d 2 P{K)/da 2 0. For comparison with the level 
clustering (dimerization) found in one-boson spin 
models ( 8 21 ) it is interesting to demonstrate 
the numerical evidence of the trimerization in the 
E (gi e model. From Fig. 8 it is evident that 
the clustering of levels with two and three dom- 
inant spacings at fixed angular momenta causes 
notable deviations from both Poisson and Wigner 
distributions. The third (localized) peak is char- 
acteristic for the two-boson case of a two-level 
model, while the one-boson case is distinguished 
by the two peak structure. The long-range statis- 
tics comprises a characteristics of the spectrum 
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p < s > a=2, j=S3/2 



o.e o.s 1 i.; n o.;5 0.5 0.75 1 1.25 1,5 i«Ts 



0.5 1 1.5 



■ 



0.5 1 1.5 




Figure 8. Level spacing distributions for partic- 
ular j. Dimerization clustering related to two 
dominating level spacings is evident for a = 2 
and j = 63/2; Third peak from localized states 
occurs at large j and a, dominating close to the 
semiclassical limit (compare with Fig. 4). The 
crossover from the linear dimerized regime to the 
nonlinear localized regime with increasing j and 
a is evident. 



defined as its " rigidity measure" . Among possible 
forms of this parameter one can name the fluctua- 
tion of level number 5 2 N in an energy band of the 
widths E, the A 3 statistics |16I17I18| etc. The 
RMT, in its turn, predicts a characteristic be- 
haviour for all these measures which should scale 
as log-E (or log (TV)), meanwhile for fully uncorre- 
lated sequences of levels the linear scaling ~ (N) 
is expected. In Fig. 8 the sample example of the 
long-range statistics (A3) is given. It is seen that 
the behaviour of the spectral rigidity of the JT 
systems conforms rather to the RMT patterns. 

4. Vibronic spectra of E ® e Jahn-Teller 
model as a Calogero-Moser gas of 
pseudo-particles with repulsion. 

The approach of the present section is to some 
extent complementary to the statistical approach 
of the preceding section. The mapping of en- 
ergy spectra E n (a) of non-integrable Hamiltoni- 
ans with one parameter of non-integrability (a), 
H = Ho + aV where Ho is integrable, on a clas- 
sical integrable Calogero-Moser model of inter- 
acting pseudo-particles moving in a pseudo-time 



Figure 9. Example of the A3-statistics (spec- 
tral rigidity) for the JT model. The long-range 
statistics appears to show the log-like behaviour 
as function of (N) typical for the RMT predic- 
tions (fully developed quantum chaos). 



a is a very useful tool for investigation of en- 
ergy spectra of non-integrable systems first in- 
troduced by Pechukas and Yukawa ^Hj- The 
mapping is based on definition of the dynamic 
variables E n (a) = x n (r), dx n /dr — p„(r), with 
the pseudo-time r = a. Respective classical 
dynamic equations for pseudo-particles evolving 
in the pseudo-time r with repulsive interactions 
can be rewritten for fluctuations 8m defined as 
X2n+i — %2n = 1 + #2n- The Calogero-Moser set 
of equations reads 



dp n 
dr 



= 2 



7, L, 



J nm JU mn 



dL r , 



dr 



where 



E 



L m lLlr_ 



,(io) 



p n = (n(T)\V\n(r)) =V nn . (11) 

V mn are matrix elements of the Hamiltonian. 
Special case is the Hamiltonian with excited spec- 
trum consisting of 2 x 2 diagonal clusters with 
diagonal and nearest neighbor off-diagonal ele- 
ments, f nrnr = = f dnr and fn r + l,n r = fon r +V 

respectively. For the case Q we have V2 n +x,2n = 
fd Ur , V 2n ,2n-i = fo„ r - There is 

L2n+I2n = fdn r ( x 2n+l (t) - IJnW), 
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f dnr = 2y/2y/n+l + \j-l/2\, 

L2n-12n = fo nr ( X 2n (t) ~ X 2n -l (t)), 



(12) 



For what follows we define small fluctuation S by 
X2n+i -x 2n = l + fan, and approximate (x 2n +i - 
x 2n )~ 1 » 1 - S 2n + - S 2n- Then, for energy 
fluctuations 5 2n , S 2n +i from J7J, (fTUf) and ltT2|) we 
get nonlinear equations with the lowest powers of 
nonlinearity in the form 



d 2 5< 



2n+l 



Ot 2 



4/d„_.+l(^2n+2 + $2r, 



25- 



2n+l 



8 (/dn r + l ~ fon r + i)S 2 n+l 



^(fdl r + l - fdi r )$2n 



8/on r +1^2n+l + 4/d nr + i^2ri+2 + ^fdn r ^2n + 



-0(£ 3 ),(13) 



4/ 2 r + 1 (5 2 Ti + l + $2n-l - 2 ^2n) 

= 8 (/o« r + l - /dn r +l)^ 2 ™ - 
4(/o^ r + l - !o n y>2n-\ 

-8fd1 lr 8l n + 4/ 0l 2 lr+1 ^2„ +1 + 4/ 2 lr ^„_ 1 

+ 0(£ 3 ). (14) 

Eqs. (|13[1 and (|14|l can be approximated by tak- 
ing second terms on the l.h.s. as second deriva- 
tives. In view of definitions (|12|) for large j 3> n 
we have \f OJkr+1 < \fd nr \, fd(n r + 1) « h rw - If 
we redefine <5 2n — S2n — 1 and assume for the de- 
viations plausible relations b\ n rs S2n+2 ~ ^In+i 
in Eq. lj*H?)l and (|H|) . we obtain, up to the second 
order 



d 2 * 



2n+l , p 2 d 2 S 2n+ i 

- 4/d nr+ r 



dr 2 



dn 2 



8fdl r+ i(52 



n+l 



J 2n+lh 



(15) 



and imaginary space, r — > it, n — > in, respec- 
tively. 

Analogously, alternative set of equations reads 



9 2 52n+l . „ 2 d 2 S 2n +i 

~~ 4 /o„ +i" 



dr 2 

d 2 S 2n 
dr 2 ' 



dn 2 



8/Z +l(^n+^„),(17) 



2 ^ 2 ^2» 

dn 2 



-8/d 2 lr + i((52n 



-^n+l)-(18) 



Transition to the imaginary space changes the 
sign of the potential, i.e. the particle turns to 
the tunneling domain and the pseudo-particles 
change their parity, e.g. hn+iivi = 2f drir -> 
$2n{vi), 5 2n {v 2 = 2f drir -> S 2n +i(v 2 = 2f 0rir+1 ), 
so that the trajectories interchange their veloc- 
ities, _S 2n (v 2 ) = 2f 0nr+1 — > 6 2n {vi) = 2f drir 
and S 2n +i(vi = 2fd„ r ^+1(^2) = 2/ „ r+1 . 
By another words, the particles interchange their 
velocities, v\,v 2 when transferring from real to 
imaginary space and vice versa. This scenario 
can be understood as series of soliton collisions 
related to two subsequent levels. 
Indeed, equation i(T5|) for large j (fd n +i ~ 
2y / \2j — 1|) admits solution in form of the prop- 
agating nonlinear pulse S(( = n — vt) 



<Wi(C ~ Co) = - -cosh 2 



C-Co 



(19) 



where (0 = no — vtq restores the translational 
invariance in the system of levels and L = (v 2 — 
(fdn r +i) 2 ) 1/2 /[2V2\fdn r +i\] is the soliton width. 
The corresponding coordinate is identified from 
dx 2n +i/dn = 1+^2,1+1 as a soliton (kink) solution 



^2n+l 



3L 

2n tanh 

2 



C-Co 



(20) 



If in equations (|15|) and i|16|l one approximates 
^2«+2 ~ 5 2n and assumes 5 2n +\ = -S 2n one ob- 
tains approximately 



dr 2 



Af 2 d 2 5 2n _ 
K+1 dn 2 ~ 



(16) 



In view of different signs of the potentials on the 
r.h.s. of equations l|15|) and ll(jfl even and odd 
pseudo-particles can be related as moving in real 



fd n , 



J2n+1 



+ 1 



fd 



2fon r 



"<52n(l - S 2n )- 



(21) 



Equation (|21|l is a well known logistic equation 
implying transition to the chaotic regime for A = 

fd \tl f + f r ~ ra+|j ~ 1/21 > A crit = 3.56994 . . .. 
Conventionally, the term "quantum chaos" is 
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used to denote the traces of classical chaotic be- 
havior at a quantum level. As already men- 
tioned in the Introduction, the classical coun- 
terpart of the system under consideration can- 
not be defined uniquely. The semiclassical ap- 
proximation in two-level systems generically leads 
to classical chaotic patterns as a result of non- 
linear coupling between two subsystems, boson 
and electron, considered respectively as classical 
and quantum; The onset of classical chaos cor- 
responds to the energies above the first diabatic 
line. It is to be emphasized that the mentioned 
chaotic behavior refers to purely quantum regime 
of medium j and n r between the weak coupling 
with j <C n r (dimerized pairs of oscillators) and 
strong coupling with j 3> n r (kink lattice) do- 
mains. Thus, this chaotic behavior can be re- 
garded as being of essentially quantum nature. 
One can conclude, that the mapping of the quan- 
tum system on the classical Calogero-Moser gas 
with repulsive interactions enables one to use the 
classical formalism for describing the system via 
its quantum numbers. So the remarkable feature 
of this approach is its ability to represent quan- 
tum chaos by classical equations. 
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